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We report extensive all-electron time-dependent density-functional calculations and nonresonant 
inelastic x-ray scattering measurements of the dynamical structure factor of 3d transition metals. 
For small wave vectors, a plasmon peak is observed which is well described by our calculations. 
At large wave vectors, both theory and experiment exhibit characteristic low-energy electron-hole 
excitations of d character which correlate with the presence of d bands below and above the Fermi 
level. Our calculations, which have been carried out in the random-phase and adiabatic local-density 
approximations, are found to be in remarkable agreement with the measured dynamical structure 
factor of Sc and Cr at energies below the semicore onset energy (M-edge) of these materials. 

PACS numbers: 78.70.Ck, 71.15.Mb, 71.45.Gm 



I. INTRODUCTION 

Electron-hole and plasmon excitations in solids have 
been investigated for many years on the basis of the 
wave- vector and frequency dependent dynamical electron 
density response of the solid. ^ The first theoretical in- 
vestigations of the dynamical density response function 
which take into account explicitly the band structure 
of the solid were performed during the last decademi^ 
Since then, first-principles calculations of the dynamical 
response of real solids have been reported for simple met- 
als and semiconductors , ^1^1^ and also for the noble met- 
als Cu)^ AgjS and Au,^" a few transition metalsfii*i^ the 
transition-metal oxides Ti02^^-^'' and NiO^ MgB2r& 
and manganitesiii Most of these calculations discuss the 
response of the solid in the framework of time-dependent 
density functional theory (TDDFT) within the random- 
phase approximation (RPA) and also incorporating the 
many-body effects, i.e. exchange and correlation effects, 
using some local-density approximation (LDA). In gen- 
eral, these approximations have proven to be quite ef- 
fective when studying excitations which do not involve 
localized states. 

Vast et al.^^ have recently reported ab initio calcula- 
tions of the electron energy-loss spectrum of the transi- 
tion metal oxide Ti02, focusing on the excitation from 
the Ti 3p semicore levels. They have compared their re- 
sults with electron energy-loss spectroscopy (EELS) mea- 
surements and have found that the inclusion of crystal 
local-field effects (CLFE), which arise when the micro- 
scopic electric field varies rapidly over the unit cellji^ 



turns out to be crucial for a faithful description of the 
experimental spectrum for small wave vectors. 

The dynamical electron density response of solids 
can be probed with the use of optical absorption 
experiments,^^ EELSjiSiSS and inelastic x-ray scatter- 
ing (IXS) techniquesiSi In particular, the advent of syn- 
chrotron sources and high-resolution IXS facilities have 
allowed detailed investigations of the processes respon- 
sible for short-range many-body effects at large wave 
vectors iSSiSiSi Macrander et alm^ have reported the first 
IXS measurements showing a strong spectral peak for 
excitations from Ti 3p semicore states to empty states 
above the Fermi level. In particular, these authors mea- 
sured the dynamical structure factor for Ti and TiC sin- 
gle crystals for intermediate and large momentum trans- 
fers, and in the case of TiC they compared their measure- 
ments to local-density-functional calculations with no in- 
clusion of crystal local-field effects. Montano et al?^ in- 
cluded these effects in the study of the dynamical struc- 
ture factor of SiC for large wave vectors, finding the in- 
clusion of crystal local-field effects to be very important 
in order to explain the measured IXS spectrum. More 
recently, Gurtubay et alm^ have found that large crys- 
tal local-field effects in the dynamical structure factor of 
Ti02 yield a sharp loss peak in the low-energy regime 
(~14 eV), which is also present in the experimental IXS 
spectra. 

In this paper, we report extensive all-electron TDDFT 
calculations and nonresonant IXS measurements of the 
dynamical structure factor of various 3d transition met- 
als with d states below and above the Fermi level. For 
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small wave vectors, a plasmon peak is observed which is 
well described by our calculations. At large wave vectors, 
both theory and experiment exhibit a sharp peak which 
is absent for small wave vectors and which we argue to 
be originated in the presence of direct d-to-d transitions 
involving d states below and above the Fermi level. Com- 
parison between theory and experiment indicates that at 
low energies the basic physics of the dynamical response 
is contained in the random-phase approximation (RPA) 
as long as the crystal local-field effects are taken into 
account. Furthermore, our results indicate that the in- 
clusion of many-body effects through an adiabatic LDA 
(ALDA) kernel yields an unexpectedly good description 
of the electronic response of valence electrons of predom- 
inant d character. 



II. THEORY 

A. Single-particle Bloch states: LAPW basis 

The starting point of our calculations is a set of well- 
converged Bloch states (j)k,n and energies ek,n, which are 
the eigenf unctions and eigenvalues of the Kohn-Sham 
equation of density-functional theory (DFT).^^'^* For an 
accurate description of the dynamics of transition metals 
with narrow d bands, we expand these Bloch states in a 
linearized augmented plane wave (LAPW) basis, ^^'^^ as 
follows 



ik,„(r)=5]Ck,„(G)V'LAPW 



(r), 



(1) 



G 



G being vectors of the reciprocal lattice. 

The choice of the basis functions V'k+G^('') made by 
partitioning the primitive unit cell into non-overlapping 
muffin-tin spheres centered on the atomic nuclei, where 
the electronic environment retains an atomic-like nature, 
and the interstitial region, where electrons may be de- 
scribed by plane waves. The full effective potential enter- 
ing the Kohn-Sham equation of DFT is expressed accord- 
ingly, i.e., it is expanded in lattice harmonics inside the 
atomic spheres and it is described by using a Fourier se- 
ries in the interstitials. The maximum reciprocal-lattice 
vector Gmax that we consider in the expansion of Eq. 
will be dictated by the cut-off parameter Rmt x Gmax, 
Rmt being the radius of the smallest muffin-tin sphere 
in the unit cell. 



B. Kohn-Sham density-response function 

Once we have an accurate description of the Kohn- 
Sham single-particle Bloch states (/)k,n and energies £k,n, 
we evaluate the density-response function x''(r,r',w) of 
noninteracting electrons moving in the effective Kohn- 
Sham potential of DFT.^^ Since this quantity must keep 



the periodicity of the solid, one can write 

BZ 



k G,G' 



(2) 

where k is a vector in the first Brillouin zone (BZ) and the 
Fourier coefficients Xg c (^j ^) given by the following 
expression: 



^ BZ 



XG,G'(k,t^) 



/>k',n|e~''^*'+'=^^-'"|'/'k'+k,n')(0k'+k,„'|e 



/k'+k,i 



k'+k,n' + HtO + il]) 

■'(•'+'''^-""|0k',„).(3) 



Here, n represents the normalization volume, 77 is a pos- 
itive infinitesimal, and /k.n are Fermi factors, which at 
zero temperature (T = 0) are simply given by the Heav- 
iside step function 



/k,n — 2 Q{eF — £k,n), 



(4) 



ep being the Fermi energy of the solid. 

The major task in the calculation of the Fourier co- 
efficients XQQ,(k, w) resides in the computation of the 
matrix elements entering Eq. Q. Using Eq. IpJ, both 
the integration over angular variables and the sum over 
reciprocal-lattice vectors can be done analytically, and 
one finds explicit expressions for the matrix elements in 
terms of one-dimensional integrals involving basis radial 
functions and spherical Bessel functionsi^ 



C. Interacting density-response function 

In the framework of TDDFT,22i2i^ the exact density- 
response function x(r, r', ^) of a periodic solid of electron 
density no(r) can be expanded in a Fourier series of the 
form of Eq. with the Fourier coefficients XG,G'(k, uj) 
given by the following expression: 

XG,G'(k,w) = XG.G'(k,t^) + XG,Gi(k,t^) 

Gi,G2 

X {i;Gi(k)(5Gi,G2 +/G?,G2M(k,w)}xG2,G'(k,w), 

(5) 

where Xg G'(k, i^) are the Fourier coefficients of Eq. 
vgO^) = 47re^/|k -I- Gp is the Fourier transform of 
the bare Coulomb potential, and /q*^, [n](k, w) are the 
Fourier coefficients of the functional derivative of the 
time-dependent exchange-correlation (XC) potential of 
TDDFT. In the RPA, the XC kernel f^%[n]{k,uj) is 
taken to be zero; in the ALDA, it is approximated by an 
adiabatic kernel of the form 



dVxcin) 



dn 



(6) 



n=no(r) 



where Vxcin-) is the XC potential of a homogeneous elec- 
tron gas of density n. 
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D. Dynamical structure factor and x-ray inelastic 
scattering cross section 

The dynamical structure factor S{q,u}) of a many- 
electron system, which determines the electron-density 
fluctuations at a given wave vector q and frequency cj, is 
connected to the imaginary part of the density-response 
function r',w) through the fluctuation-dissipation 
theoremi^ In the case of a periodic solid at T = 0, one 
finds 

S{k + G,u;) = -2hmmxGMKi^), (7) 

k being a wave vector in the first BZ. 

Within the first-Born approximation, the double differ- 
ential scattering cross-section for x rays to transfer mo- 
mentum ?i(k + G) and energy huj to a periodic solid is 
given by 

where (e^, e/) and (wi, tof) refer to the polarization vec- 
tor and frequency of the incident and scattered photon, 
respectively. 

E. Dielectric matrix and macroscopic dielectric 
function 

In order to investigate the contribution to the energy- 
loss spectrum coming from the excitation of plasmons, it 
is convenient to consider the wave-vector and frequency 
dependent inverse dielectric matrix of a periodic solid 
which yields the potential induced at a given point of 
space by an arbitrary (small) frequency-dependent ex- 
ternal potential. One finds 

'^G',G'i^^^) = ^G,G' +UG(k)xG,G'(k,w), (9) 

with k being a vector in the first Brillouin zone. 

For k wave vectors in the first Brillouin zone, the 
macroscopic dielectric function is defined as follows 

eAf(k,c^) = lAo,o(k,c^), (10) 

where e(^J(k, cj) represents the first diagonal element of 
the inverse dielectric matrix of Eq. 0. In the absence 
of crystal local-field effects (CLFE), the off-diagonal ele- 
ments of the inverse dielectric matrix are negligible and 
the macroscopic dielectric function coincides, therefore, 
with the first diagonal element of the dielectric matrix, 
eG,G(k, (x)), which in the RPA is given by the following 
expression; 

eg^c^,(k,c.) = 5g.g' -«G(k)x^,G'(k,c.). (11) 

Collective excitations (plasmons) are known to occur 
at energies where both the real and the imaginary part of 



the macroscopic dielectric function are close to zero. On 
the other hand, the diagonal elements of the imaginary 
part of the RPA dielectric matrix are known to repre- 
sent a measure of the number of states available for real 
transitions involving a given momentum transfer ?iq and 
energy transfer huj. 

III. EXPERIMENT 

Nonresonant inelastic x-ray scattering measurements 
of the dynamical structure factor ^(q, lo) of single crys- 
tals Sc and Cr have been obtained using the UNICAT 
undulator Beam Line on the Advanced Photon Source 
(APS) at Argonne National Laboratory. The measure- 
ments were made in symmetric refiection geometry, with 
7.6 keV x-rays, a range of wave vectors from 0.5 to 3.98 
A , and an energy resolution of 1.1 eV. The shape- 
asymmetry of the tails of the quasi-elastic scattering near 
zero energy transfer for the analyzer and scattering ge- 
ometry used was determined directly by measurements 
made on wide band-gap (> 10 eV) insulators (e.g. LiF, 
CaF2), which indicated the tails to be symmetric for 
losses greater than 2.5 eV. The quasi-elastic contribution 
to the scattering for Sc and Cr was then removed using 
measurements on the energy gain side of the quasi-elastic 
peak and restricting the data to energy losses of 2.5 eV 
or greater. The rapid falloff of the quasi-elastic contri- 
bution rendered the tails insignificant for energy losses 
greater than 6-8 eV for the conditions used. 

In order to reduce to absolute units the APS spectra 
of the transition metals Sc and Cr, we have first per- 
formed an absolute measurement on Al.^^ For measure- 
ments that do not extend past the Al L-edge at ~ 72 eV, 
only the three valence electrons of this metal contribute 
to the first-moment / sum-rule 

/ dww s(q,w) = a (12) 
Jo "^e 

Here, uq is the average electron density of the 3s^ 
valence electrons in Al, and s(q, w) is the reduced dy- 
namical structure factor: s(q, w) = S {q, oj) / {hfl) . Ac- 
cording to Eq. JHl) , the dynamical structure factor should 
be proportional to the relative intensity /(q, a;)//o of the 
scattered beam, Iq being the incident beam power; hence, 

s(q,c^) = C/(q,c^)//o. (13) 

The maximum energy Tiujmax entering the left-hand side 
of Eq. H12|l must be below the 72 eV semicore onset en- 
ergy and the coefficient a represents the relative con- 
tribution of valence 3s^ ip^ electrons to the / sum-rule 
for energies below humax- We have carried out ab ini- 
tio pseudopotential calculations of the dynamical struc- 
ture factor of Al, and for htOmax = 65 eV we have found 
a = 0.95. 

Introducing Eq. (|13|) into Eq. H12|) one finds the cal- 
ibration factor C, which we have derived from our IXS 
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measurements of the I{q,uj)/Io ratio in an Al sample. 
(See Ref. Is^ for details.) Since this calibration factor is 
inversely proportional to the effective volume of material 
sampled by the x-ray beam (and is independent of the 
atomic configuration of the material) and for thick sam- 
ples and symmetric Bragg geometry the volume ratio of 
the material sampled in different targets is given by the 
ratio of the corresponding linear absorption coefRcients, 
the dynamical structure factor sx (q, w) of an arbitrary 
(X) target can be obtained from the following expression: 

sx(q, w) = C/ix/ma; /x(q,t^)//o, (14) 

with fix and /iA; being the linear absorption coefficients 
of the X material and Al, respectively. 

As suggested in Rcf. 36, analogous (albeit sample ge- 
ometry and q-dependent) expressions for the ratio of the 
effective scattering volumes can be derived for transmis- 
sion geometry and non-symmetric reflection scattering 
geometry. 

IV. RESULTS AND DISCUSSION 

In this section, we first present TDDFT calculations 
and nonresonant IXS measurements of the dynamical 
structure factor of the hexagonal close-packed (hep) Sc 
and the body-centered cubic (bcc) Cr. Then, we present 
a comparison of the dynamical response of the transition 
metals Sc, Ti, V, Cr, Fe, Co, and Ni, all belonging to the 
3d series, and we focus on the evolution of a low-energy 
peak that is visible at large wave vectors and which is 
originated in the presence of d-to-d transitions involving 
d-states below and above the Fermi level. 

All the calculations presented below have been car- 
ried out by first expanding Bloch states in an LAPW 
basis, as in Eq. with a cutoff parameter of Rmt x 
Gmax — 8, using spherical harmonics inside the atomic 
spheres up to / = 10, and including 3s and 3p states as 
semicore states. We have then solved self-consistently 
the Kohn-Sham equation of DFT in the local-density 
approximation (LDA), with use of the Perdew-Wang 
parametrizationSi of the Ceperley-Alder XC energy of 
a uniform electron gas^^ we have calculated the Fourier 
coefficients Xg g 0^^ ^) from Eq. Q by keeping unoccu- 
pied KS states up to 7.5 Ry, and we have finally solved 
Eq. lO either in the RPA or the ALDA with a finite 
damping parameter 77 = 0.65 eV. 

In the case of the hexagonal close-packed Sc, Ti, and 
Co, the coefficients XG,cj(k, w) have been evaluated for 
wave vectors perpendicular to the hexagonal plane, i.e., 
along the (001) direction, and well-converged results have 
been obtained by using a 8 x 8 x 16 sampling of the BZ 
corresponding to 90 points in the irreducible BZ (IBZ). 
In the case of the body-centered cubic Cr, V, and Fe, cal- 
culations have been performed for wave vectors along the 
(100), (110), and (111) directions, and well-converged re- 
sults have been obtained by using a 10 x 10 x 10 sampling 
of the BZ corresponding to 47 points in the IBZ. Finally, 



for the face-centered cubic Ni the coefficients XG,G(k, lu) 
have been evaluated for wave vectors along the (100) di- 
rection, and well-converged results have been obtained 
by using a 12 x 12 x 12 sampling of the BZ corresponding 
to 72 points in the IBZ. 

In order to investigate the impact of crystal local-field 
efects, we have first neglected these effects by considering 
only the diagonal elements of the Kohn-Sham density- 
response matrix entering Eq. ((SJ (diagonal calculation), 
and we have then solved the matrix equation with a given 
number of G vectors (full calculation). Well-converged 
results have been obtained with the use of 35, 43 and 51 
G vectors in hep, bcc and fee structures, respectivelyi^^ 

A. Sc 




energy (eV) 

FIG. 1: (Color online) Density of states (DOS) of Sc, at en- 
ergies near the Fermi level (dotted line) for each component 
of the angular momentum inside the atomic sphere. Dash- 
dotted black, dashed red and solid green lines denote the 
number of states per atom and per eV with s, p, and d char- 
acter, respectively. The inset shows a wider energy range, in 
which the 3s (dash-dotted black line) and 3p (dashed red line) 
semicore states can be distinguished. 

Sc ([Ar] 3d^ 4s^) being the first of the transition metals 
has a single 3d electron per atom. Hence, the 3d band 
is located below and above the Fermi level, as shown by 
the density of states (DOS) drawn by a solid green curve 
in Fig. n The inset of this figure shows the location of 
the 3s^ (dash-dotted black line) and 3p^ (dashed red line) 
semicore bands, which are located at ~49 eV and ~28 eV 
below de Fermi level, respectively. 

Early measurements of the energy- loss spectra of high- 
energy electrons in Sc were reported by Brousseau- 
Lahaye et al. for small wave vectors?^ A plasmon peak 
was identified at ~ 13 eV, which agrees with the ex- 
pected plasmon energy of valence electrons (3d^4s^) in 
Sc. Pseudopotential-based first-principles calculations 
of the energy-loss function of Sc were first reported by 
Schone and Ekardliii at small and intermediate wave vec- 
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tors and more recently by Gurtubay et al.i^ A detailed 
analysis of the plasmon energy dispersion of this material 
was reported in Ref. |3 demonstrating that the plasmon 
peak at ~ 13 eV corresponds indeed to zero values of 
the real part of the macroscopic dielectric function where 
the imaginary part is small and showing that the crystal 
local-field effects have a considerable impact on the plas- 
mon energy at small wave vectors which gives rise to an 
interplay with the XC effects built into the many-body 
kernel. 
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FIG. 2: (Color online) Reduced dynamical structure factor 
s(q, u) = ^(q, uj) I (MX) of Sc, as obtained in the RPA (dashed 
black line), ALDA (solid red line), and by IXS measurements 
(green circles). The wave vector has been taken to be q= 
8/16(001) (left panel) and q= 12/16(001) (right panel) in 
units of 27r/c, c being the lattice constant in the 2-direction. 

In Fig. 12 we display new first-principles calculations 
(dashed black and solid red lines) and IXS measurements 
(green circles) of the dynamical structure factor s(q, 

of Sc at wave vectors of mag nitude |q|=0.6 A (left 

panel) and |q| = 0.89 A (right panel) and for energies 
below the ~ 28 eV semicore onset energy (M-edge) of this 
material. At these wave vectors, the dynamical structure 
factor displays a broad plasmon peak at ~ 15 eV, which 
corresponds roughly to the expected plasmon energy of 
three valence electrons per atom in Sc (r^ — 2.38);'*^ nev- 
ertheless, the shape of this plasmon peak, which mixes 
with the fine structure caused by d-electron transitions, 
bears little resemblance to the plasmon peak of a ho- 
mogeneous electron gas. In particular, some structure 
is present at energies between 5 and 10 eV; this struc- 
ture, which is also present in the imaginary part of the 
dielectric matrix and has also been predicted to exist for 
smaller wave vectors, ^^'^^ occurs in a region where inter- 
band transitions from the valence bands to unoccupied p 
states above the Fermi level take place. 

Fig. [3 also shows that the agreement between our ah 
initio calculations and the experimental measurements is 
excellent, both in the energy of the plasmon and in the 
overall shape of the energy-loss spectrum. Furthermore, 
for energies below the plasmon peak the agreement be- 
tween ALDA calculations (solid red line) and IXS mea- 
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FIG. 3: (Color online) As in Fig. H but for q= 12/16(001) 
(in units of 2ti/c) and energies that include the excitation of 
3p semicore states. The blue dotted line denotes the diagonal 
(i.e. without CLFE) RPA calculation. 



First-principles calculations (dashed black and solid 
red lines) and IXS measurements (green circles) of the 

dynamical structure factor of Sc at |q|=0.89 A ^ are 
exhibited in Fig. |21 for energies up to 70 eV, well above 
the ~ 28 eV semicore onset energy (M-edge). For com- 
parison, diagonal RPA calculations are also plotted in the 
same figure (blue dotted line), showing that the impact 
of crystal local-field effects is very large just above the 
M-edge, where transitions from the occupied Sc semicore 
3p states to the lowest conduction bands occur. Crystal 
local-field effects bring the onset of the M-edge to larger 
energies and reduce the corresponding peak height, in 
close agreement with experiment. The shift to larger 
energies stems from the fact that the semicore states in- 
volved in these transitions are highly localized in real 
space and therefore neglecting the off-diagonal terms in 
Eq. and subsequently in Eq. @, underestimates the 
large dynamical screening effects near the semicore edge. 

Figure O shows, however, that although our full RPA 
calculation (black dashed line) is in very good agreement 
with experiment (green circles) there is a mismatch be- 
tween the ALDA calculation (red solid line) and the mea- 
sured spectrum. This discrepancy calls for further inves- 
tigations of (i) an improved description of single-particle 
energies beyond the LDA and (ii) accurate many-body 
kernels accounting for XC effects beyond the ALDA. 
Both RPA and ALDA calculations have been carried out 
with the very same form of the XC potential entering 
the KS equation of DFT. Hence, we can rule out our 
approximate XC potential as the origin of the large dis- 
crepancy between ALDA calculations and experiment, 
and sources of error of our ALDA calculations should 
be searched both in the adiabaticity and the locality of 
the ALDA XC kernel in a solid where inhomogeneities 
of both the static unperturbed electron density and the 
induced density play an important role, especially when 
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FIG. 4: (Color online) Dynamical structure factor s(q,tj) of 
Sc, as obtained in the RPA (black dashed line), ALDA (red 
solid line), and by IXS measurements (green circles). The 
wave vector has been taken to be q= [4/16(001) + (002)] (up- 
per panel) and q= [8/16(001) + (002)] (lower panel) in units 
of 2-k/c, c being the lattice constant in the z-direction. 



semicore states with large binding energies are involved. 

In Fig. 01 we show calculations and measurements of 
the dynamical structure factor of Sc, as in Fig. 2, but 
now at wave vectors of magnitude |q| = 2.09 A (upper 

panel) and |q|=2.98 A (lower panel). At these large 
wave vectors, the plasmon energy falls into the electron- 
hole pair excitation continuum so there is no well-defined 
plasmon peak. Instead, a new feature appears in the dy- 
namical structure factor at ~ 4eV, which is very sensi- 
tive to crystal local-field effects. A close examination of 
the various transitions that contribute to the KS density- 
response function of Eq. @ indicates that this new fea- 
ture, which is absent for small wave vectors, is originated 
in the presence of d-to-d transitions involving narrow d 
complexes lying below and above the Fermi level. Such 
direct d-to-d excitations have not been observed in opti- 
cal absorption measurements (where the wave vector is 
zero) and EELS experiments (where the wave vector is 
small), which underscores the importance of IXS for the 



examination of the electronic structure of these materi- 
als. 

An inspection of Fig. ^ also shows that introducing 
many-body XC effects within the ALDA (red solid line) 
leads to a faithful description of the low-energy part of 
the measured energy- loss spectrum (green circles), as oc- 
curs for smaller wave vectors (see Figs. |2| and IS)), even 
at ~ 4 eV where the energy loss is dictated by the pres- 
ence of d-to-d transitions. At energies above the M-edge, 
the agreement between our full RPA calculations (black 
dashed line) and the IXS measurements is excellent, but 
there is a considerable mismatch between the ALDA cal- 
culation (red solid line) and the measured spectrum, as 
occurs for smaller wave vectors (see Fig.l^J. 

A closer look into the semicore energy region in Figs.|3| 
and 0] still shows a 2-3 eV mismatch between the cal- 
culated RPA and the measured IXS semicore excitation 
onset. This stems from a well-known deficiency of the 
local-density approximation we have used to describe our 
single-particle orbitals^^ which in general underbinds the 
energy of highly localized occupied states and is known 
to occur even in the case of a simple metal such as Ali^ 



Cr 



Cr ([Ar] 3d^ 4s^) lying in the middle of the 3d period 
exhibits a high density of both occupied and unoccupied 
d states near the Fermi level. Hence, this transition metal 
is the ideal candidate for the investigation of the char- 
acteristic d-to-d transitions that we have observed in the 
case of Sc. 
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FIG. 5: (Color online) Density of states (DOS) of Cr near the 
Fermi level (dotted line) for each component of the angular 
momentum inside the atomic sphere. Black dash-dotted, red 
dashed, and green solid lines represent the number of states 
per atom and per eV of s, p, and d character, respectively. 
The inset shows a wider energy range, in which the 3s (black 
dash-dotted line) and 3p (red dashed line) semicore states are 
visible. 

The density of states of Cr is depicted in Fig.|31 where 
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d states distributed at both sides of the Fermi level are 
represented by green solid lines. As in Fig. ^ the inset 
of this figure shows the location of the 3s (black dash- 
dotted line) and 3p (red dashed line) semicore states at 
~70 eV and ~42 eV below the Fermi level, which hold 2 
and 6 electrons per atom, respectively. 

Inelastic x-ray scattering measurements of the energy- 
loss spectra of Cr were reported recently by Montano and 
Macrander^^ for momentum transfers ranging from 0.80 
to 5.0A along the (110) direction. Spectral features 
were found corresponding to plasmon excitations, as well 
as a second peak corresponding to excitations from the 
3p semicore states to the 3d valence states; however, no 
structure was visible at low energies near ^ 4 eV, where 
the x-ray spectrum was dominated by the quasielastic 
phonon scattering feature. To the best of our knowledge, 
no theoretical investigations of the dynamical response 
of this material have been reported before. 




; 
; 
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FIG. 6: Top panel: RPA full (solid line) and diagonal (dotted 
line) calculations of the dynamical structure factor s(q,tj) of 
Cr. Bottom panel: RPA calculations of the real (dashed line) 
and imaginary (solid line) parts of the first diagonal element 
of the dielectric matrix eG,G'(q,'^) of Cr. In both panels, the 
wave vector has been taken to be q= 1/10(110) in units of 
2n/a, a being the lattice parameter. 

The top panel of Fig.|Blshows first-principles RPA diag- 
onal (dotted line) and full (solid line) calculations of the 
dynamical structure factor s(q, w) of Cr for small wave 
vectors of magnitude |q| = 0.15 A~^ along the (110) di- 
rection. At energies below the ~ 43 eV semicore onset en- 
ergy (M-edge), where crystal local-field effects are found 
to be small (diagonal and full calculations are very close 
from each other), we find a prominent broad feature at 



26 eV, which agrees with the expected plasmon energy 
of 6 valence electrons per atom in Cr (r^ = 1.48) and is 
clearly identifiable as a collective excitation of hybridized 
3d^ and 4s ^ valence electrons. Indeed, this peak corre- 
sponds (see the bottom panel of Fig. O to a zero value 
of the real part of the matrix element eo,o(cii^^) where 
Imeo.o(q, ti^) is small.^^ 

At energies above the M-edge, transitions from the lo- 
calized Cr semicore 3p states below the Fermi level to 
the lowest conduction bands take place. As a result, the 
impact of crystal local-field effects is large (see the top 
panel of Fig. EJ, as occurs in the case of Sc. These ef- 
fects bring the onset of the M-edge to larger energies and 
considerably reduce the peak height above the M-edge, 
which becomes a broad feature around 50-60 eV involving 
primarily excitations out of the 3p band of Cr. 




Q)(eV) 

FIG. 7: (Color online) RPA calculations of the dynamical 
structure factor s(q,tj) of Cr at wave vectors of magnitude 
|q| ~ 0.8 along three high-symmetry directions. Dotted 
green, dashed red, and solid black lines correspond to q= 
4/11(100), q= 5/20(110), and q= 2/10(111) wave vectors, 
respectively, in units of 2-71 /a, a being the lattice parameter. 
The inset shows the real (dashed line) and imaginary (solid 
line) parts of the first diagonal element of the RPA dielectric 
matrix ecc (q,6j) at |q| = 0.75A . 

In order to investigate the dependence of the dynam- 
ical structure factor of Cr on the direction of the wave 
vector, we have plotted in Fig. [7| the RPA dynamical 
structure factor s(q, w) of Cr at wave vectors of mag- 
nitude |q| « 0.8 A-i along the (100), (110), and (111) 
directions^! This figure shows that the dynamical struc- 
ture factor of Cr is rather similar for all high-symmetry 
propagation directions. In particular, we note that the 
long-wavelength collective excitation at ^ 26 eV (see 
Fig. 0) is still visible along all propagation directions 
when the magnitude of the wave vector is |q| « 0.8 A~^. 
This is in agreement with the IXS measurements re- 
ported by Montano and Macrander— for this magnitude 
of the wave vector along the (110) direction. Further- 
more, that this feature is still originated in the build- 
ing up of a collective excitation follows from the inset of 
Fig. where the real and imaginary parts of the ma- 
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trix element eo,o(q, t^) are displayed: this feature corre- 
sponds to a zero value of Reeo.o(q, ^) (dashed line) where 
Imeo,o(q, f^) (solid line) is small. On the other hand, 
a comparison of the high-energy part of the energy-loss 
spectra of Figs. and [7| indicates that the high-energy 
broad feature around 50-60 eV, which is also rather in- 
sensitive to the propagation direction and which involves 
primarily excitations out of the ip band of Cr, has a con- 
siderable wave-vector dependence, becoming more pro- 
nounced as the wave vector increases. 
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FIG. 8: (Color online) Dynamical structure factor of Cr 
(CLFE, i.e. off-diagonal terms in Eq.lO, included) along the 
(111) direction for wave vectors q — 6/11(001) (top panel) and 
q = 10/11(001) (bottom panel), within RPA (black dashed 
line) and ALDA (red solid line). The blue dotted line rep- 
resents the calculation without CLFE and green circles, IXS 
measurements. 

First-principles calculations and IXS measurements of 
the dynamical structure factor s(q, w) of Cr are dis- 
played in Fig. IHI at large wave vectors of magnitude 
|q| = 2.06A~^ and |q| = 3.5A~^ along the (111) di- 
rection. As in the case of Sc, at these large wave vectors 
the prominent broad peak of Figs. El and at ~ 26 eV, 
which is in the nature of collective excitations, is now 
replaced by a less pronounced broad feature at energies 
below the M-edge originated from single electron-hole ex- 



citations, and the high-energy feature above the M-edge 
is found to gain weight as the wave vector increases. 

For small wave vectors, crystal local-field effects only 
play a role at energies above the M-edge (see Figs. El 
andO, where the energy- loss spectrum is dominated by 
transitions from very localized 3p bands below the Fermi 
level. The striking feature of Fig. (HI is precisely the fact 
that as the wave vector increases huge crystal local-field 
effects are present not only at energies above the M-edge 
but also at very low energies below 10 eV. In the absence 
of crystal local- field effects (blue dotted Hue), the energy 
loss spectra at the large wave vectors under consideration 
are clearly dominated by two distinct features at 4-6 eV 
and 40-60 eV. A close examination of the various transi- 
tions that contribute to the KS density-response function 
of Eq. ^ indicates that these features involve excitations 
to unoccupied d states above the Fermi level coming pri- 
marily from occupied d and p states, respectively. As in 
the case of Sc, high-energy transitions from occupied p 
states occur for arbitrary values of the wave vector; how- 
ever, direct d-to-d transitions are not observed for small 
wave vectors. In the presence of crystal local-field ef- 
fects (dashed black and solid red lines), both peak heights 
are considerably reduced and our calculated low-energy 
feature is brought into excellent agreement with experi- 
ment. Both theory and experiment exhibit the presence 
of a very well defined low-energy peak at ~ 4 eV, which 
is originated in the presence of direct d-to-d transitions 
and had never been observed before. 

An inspection of Fig. jSl clearly indicates that (i) first- 
principles RPA calculations provide a reasonably good 
description of the measured low-energy peak at ^ 4 eV, 
provided that crystal local-field effects are included, (ii) 
inclusion of XC effects within the ALDA brings the low- 
energy peak at 4 eV into excellent agreement with 
experiment, and (iii) the ALDA, which accurately ac- 
counts for XC effects at low energies even at large wave 
vectors, fails to give a quantitative description of the 
measured energy-loss spectra at energies above the M- 
edge. This suggests that for an accurate understanding of 
the energy-loss spectra of transition metals in the whole 
energy-range under study one should go beyond and adi- 
abatic local-density approximation in the description of 
XC effects not present in the RPA. We also note that as 
in the case of Sc (see Figs. O and |3|) both RPA (black 
dashed line) and ALDA (red solid line) calculations pre- 
dict an energy onset for the excitation of 3p electrons in 
Cr that is too large, which should be a consequence of 
the failure of the LDA to describe the bonding energies 
of the semicore KS states for this material. 

The d-to-d low-energy peak that is visible at ~ 4 eV 
in Fig. El is absent for small wave vectors and is therefore 
expected to have a very strong wave-vector dependence. 
Hence, in Fig. El we have illustrated this wave- vector de- 
pendence by plotting our full ALDA calculations (lines) 
and IXS measurements (circles) of the dynamical struc- 
ture factor s(q, cj) of Cr at wave vectors of magnitude 

ranging from ^ 0.7 to S.SA along the (100) and (111) 
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FIG. 9: (Color online) Formation of the d — > d transition 
peak along the (100) and (111) directions, for different values 
of the wave vector (shown in A ). 



directions (left and right panels, respectively). First, we 
note that both theory and experiment indicate that the 
weight of the d-to-d low-energy peak at 4 eV, which 
is absent for small wave vectors, considerably increases 
with q. Secondly, it is clear from this figure that ALDA 
calculations of this low-energy feature are in remarkable 
agreement with experiment for all wave vectors and prove 



not only the presence of considerable crystal local-field ef- 
fects below the M-edge but also the importance of IXS 
for the examination of the electronic structure of a 3d 
transition metal such as Cr. 

Finally, we note that the d-to-d low-energy peak that is 
visible at ~ 4 eV in Figs. 01 and |H1 is considerably sharper 
and more intense in the case of Cr (Fig. (HJ than in the 
case of Sc (Fig. 0)). This is simply due to the fact that 
in the case of Cr the density of d states is considerably 
large not only above the Fermi level but also below, thus 
increasing the number of available d-to-d transitions. 



C. The 3d series 

As we move along the 3d transition-metal period, from 
Sc to Ni, the 3d shell fills up causing a shift of the 3d 
bands towards the occupied side of the band structure 
below the Fermi level, as illustrated by the density of 
states that have been plotted in Fig. ^| Here we focus 
on the impact that this filling up of the 3d bands has on 
the collective excitations that are present for small wave 
vectors and, more importantly, on the new low-energy 
feature that is originated in the building up of d-to-d 
transitions. 
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0.41 


0.43 


0.44 


0.37 


0.29 


large q (A ) 


2.98 


3.01 


2.91 


3.05 


3.06 
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TABLE I: Parameters used for the comparative study of the 3d transition metals. Low q and large q refer to the magnitude 
of the wave vector used in Figs. Illl and ll2l respectively. 



1. Collective excitations 

Full ALDA calculations of the dynamical structure fac- 
tor s(q, oj) of 3d transition metals are depicted in Fig. 1111 

at a low wave vector of magnitude ~ 0.4 A along the 
(001) direction and for energies below the M-edge. This 
figure clearly shows that in the case of 3c? metals with 
five or fewer d states occupied per atom (Sc, Ti, V, and 
Cr) a broad peak can be identified, which mixes with 
the fine structure due to d-electron transitions. As the 
3d band fills up (Fe, Co, and Ni), however, the number 
of states available for real transitions from occupied d 
bands below to unoccupied p states above the Fermi level 
considerably increases. Furthermore, the combination of 



band-structure effects and the building-up of collective 
modes of d electrons yields broader structures in these 
materials which may still be in the nature of collective 
excitations but which are not well defined any more. 

We note in passing that our first-principles calculations 
are in good agreement with the IXS measurements of the 
dynamical structure factor of Ti reported by Macrander 
et al.^ 



2. d d transitions 

Figure^lshows full ALDA calculations of the dynami- 
cal structure factor s(q, uj) of the 3d transition metals Sc, 
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FIG. 10: (Color online) d content of the density of states in 
the 3d transition metal row, as we move from So to Ni (in the 
figure from top to bottom and left to right). 
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FIG. 11: Dynamical structure factor of the 3d transition 
metals along the (001) direction for a similar magnitude of 
the momentum transfer (|q| ~ 0.4 A ^) and energies below 
the M-edge. The exact value of the momentum transfer can 
be found in Table Q 



V. SUMMARY AND CONCLUSIONS 



Ti, V, Cr, Fe, Co, and Ni and the corresponding IXS mea- 
surements for the 3d transition metals Sc (green circles) 
and Cr (blue squares) at a large wave vector of magni- 
tude |q| ^3 A and low energies well below the M-edge. 
This figure clearly shows the presence of a low-energy 
peak at ~ 4 eV, which follows the trends of the density 
of states depicted in Fig. ^| as should be expected from 
a feature that is originated in the building up of d-to-d 
transitions involving narrow d complexes lying below and 
above the Fermi level. The weight of this characteristic 
peak is small in Sc, which holds one single occupied d 
state per atom, it considerably increases in the case of 
Ti and V, as more d states become occupied, and it is 
extremely well defined in the case of Cr, with a half- filled 
3d band below the Fermi level. As still more d states be- 
come occupied, as occurs in the case of Fe and Co, the 
number of states available for real transitions involving 
d-to-d transitions decreases again, thereby decreasing the 
weight of the ^ 4 eV peak. Finally, in the case of Ni, with 
only two unoccupied 3d states per atom, this feature is 
hardly visible. 



We have presented a combined theoretical and exper- 
imental investigation of the dynamical structure factor 
of various 3d transition metals with d states below and 
above the Fermi level. Calculations have been carried out 
for Sc, Ti, V, Cr, Fe, Co, and Ni by first expanding single- 
particle Bloch states in an LAPW basis and then com- 
puting the density-response function of the solid either 
in the RPA or the ALDA with full inclusion of crystal 
local-field effects. The experimental dynamical structure 
factor has been obtained for Sc and Cr by first performing 
nonrcsonant IXS measurements of the relative intensity 
of the scattered beam at a given wave vector and energy 
and then reducing the measured spectra to absolute units 
from sum-rule considerations. 

For small wave vectors, a collective excitation of hy- 
bridized 3d and 4s electrons is observed, which is very 
well described by our ALDA calculations and which in 
the case of 3d metals holding five or fewer d states per 
atom (Sc, Ti, V, and Cr) has a rather well defined energy. 
At large wave vectors, IXS measurements exhibit a new 
feature at ~ 4 eV, which (i) had never been observed be- 
fore, (ii) we have found to be originated in the presence 
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FIG. 12: (Color online) Low energy peak arising from d — > 
d transitions for a large wave vector, as we move from left 
(Sc) to right (Ni) along the 3d row. The magnitude of the 
momentum transfer for each metal can be found in Table |T] 
Lines denote the ab initio ALDA calculations and circles and 
squares correspond to the normalized IXS measurements on 
Sc and Cr, respectively. 

of d-to-d transitions involving d states below and above 
the Fermi level, and (iii) is remarkably well described by 
our ALDA calculations, as long as crystal local-field ef- 
fects are fully included. This new feature is absent at the 
small wave vectors accessible by EELS, and underscores 
the importance of IXS for investigations of the electronic 
structure and the dynamical response of transition met- 
als. 

Our ALDA calculations and IXS measurements of the 
dynamical structure factor of the transition metals Sc 
and Cr have been found to show a remarkable agree- 
ment at energies below the M-edge. However, at energies 



above the M-edge significant discrepancies are found to 
exist between theory and experiment. Discrepancies be- 
tween theory and experiment at the onset of the semicore 
excitations in Cr provide a signature of the single-particle 
LDA energies underestimating the binding energy of the 
semicore states. Discrepancies at the intensity and over- 
all shape of the edge, especially in the case of Sc, should 
be due to the fact that the adiabatic local-density ap- 
proximation for the many-body kernel breaks down at 
these large energies. These discrepancies call for further 
investigations of the accuracy of various approximations 
for the XC kernel of TDDFT. 
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The density-response function of a many-electron system 
yields the electron density induced at a given point of space 
by an arbitrary (small) frequency-dependent external po- 
tential. 

The actual number of G vectors needed for convergence 
depends on the spectral feature. While the semicore is the 
feature which needs the largest numbers, fewer G vectors 
are required in the valence region. 

The electron density no of a homogeneous electron gas 
is typically characterized by the density parameter = 
(3/47rno)^^^/ao, lo being the Bohr radius. 
This deficiency contributes to the discrepancy between 
ALDA calculations and experiment; however, the large 
discrepancies that are visible in Figs. |3 and 21 should be 
searched in the adiabaticity and the locality of the ALDA 
XC kernel. 

We have carried out RPA calculations and IXS measure- 
ments of the dynamical structure factor of Al, and we have 
found that the LDA predicts the 2p semicore states of this 
material to be located at ~65 eV while the IXS measure- 
ments show the 2p semicore excitation onset at ~72 eV. 
^® At the wave vectors and energies under consideration crys- 
tal local-field effects in Cr are small and one finds, there- 
fore, eA^(q,w) « eo,o(q, c^). 

In order to be able to choose similar values of the magni- 
tude of the momentum transfer along the three directions 
under study, we have used two different samplings of the 
BZ, 10 X 10 X 10 and 11 x 11 x 11. As the damping param- 
eter rj has been taken to be = 0.65 eV in both cases, 
small differences among the three curves may be due not 
only to the fact that the magnitude of the wave vector is 
not exactly the same but also to the fact that different 
samplings of the BZ have been used with the same value 
of rj. 



